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ABSTRACT 

The atomic CII fine-structure line is one of the brightest lines in a typical star-forming galaxy 
spectrum with a luminosity ~ 0.1% to 1% of the bolometric luminosity. It is potentially a reliable 
tracer of the dense gas distribution at high redshifts and could provide an additional probe to the era 
of reionization. By taking into account of the spontaneous, stimulated and collisional emission of the 
CII line, we calculate the spin temperature and the mean intensity as a function of the redshift. When 
averaged over a cosmo logically large volume, we find that the CII emission from ionized carbon in 
individual galaxies is larger than the signal generated by carbon in the intergalactic medium (IGM). 
Assuming that the CII luminosity is proportional to the carbon mass in dark matter halos, we also 
compute the power spectrum of the CII line intensity at various redshifts. In order to avoid the 
contamination from CO rotational lines at low redshift when targeting a CII survey at high redshifts, 
we propose the cross-correlation of CII and 21-cm line emission from high redshifts. To explore the 
detectability of the CII signal from reionization, we also evaluate the expected errors on the CII 
power spectrum and CII-21 cm cross power spectrum based on the design of the future milimeter 
surveys. We note that the CII-21 cm cross power spectrum contains interesting features that captures 
physics during reionization, including the ionized bubble sizes and the mean ionization fraction, which 
are challenging to measure from 21-cm data alone. We propose an instrumental concept for the 
reionization CII experiment targeting the frequency range of ~ 200 to 300 GHz with 1, 3 and 10 
meter apertures and a bolometric spectrometer array with 64 independent spectral pixels with about 
20,000 bolometers. 

Subject headings: cosmology: theory — diffuse radiation — intergalactic medium — large scale struc- 
ture of universe 



1. INTRODUCTION 

Carbon is one of the most abundant elements in the 
Universe and it is singly ionized (CII) at 11.26 eV, 
an ionization energy that is less than that of the hy- 
drogen. With a splitting of the fine-structure level at 
91 K CII is easily excited resulting in a line emission 
at 157.7 jam through the 2 P^/2 ~^ 2 P1/2 transition. It 
is now well established that this line provides a major 
cooling mechanism for the neutral interstellar medium 
(ISM ) dDalgarno fc McCrav 1972; Tielens & H ollenbachl 
[19851 : iWolfire et all 119951 : iLehner et al.l l200l . It is 
present in multiple phases of the ISM in the Galaxy 
Wright et al. (1991) including the most diffuse regions 
Boc k et al.l (jl993h and the line emission has been de- 
tected from the phot o dissociation regions (PDRs) of 
star- forming galaxies ( Boselli et al. 2002; De Lo oze et al.l 
120111 : iNagamine et al.l 120061 : iStacev et alJ l2Q10h and, in 
some cases, even in z > 6 Sloan quasars (Wa lter et al.l 
l2QQ9h . 

The CII line is generally the brightest emission line 
in star-forming galaxy spectra and contributes to about 
0.1% to 1% of the total far- infrared (F IR) luminosity 
([Crawford et al. Ill 98 51 : IStacev et al.lll991f ). Since carbon 
is naturally produced in stars, CII emission is then ex- 
pected to be a good tracer of the gas distribution in 
galaxies. Even if the angular resolution to resolve the 
CII emission from individual galaxies is not available, 



the brightness variations of the CII line intensity can 
be used to map t he underlying distribution of galaxies 
and dark matter (|Basu et al.ll2004l : iVisbal fc Loebll201fl : 
iGong et al,ll20Tlh - 

Here we propose CII intensity mapping as an alterna- 
tive avenue to probe the era of reionization, including 
the transition from primordial galaxies with PopIII stars 
alone to star- formation in the second-generation of galax- 
ies with an ISM polluted by metals. CII intensity map- 
ping complements attempts to study reionization with 
low- frequency radio experiments that probes the 21-cm 
spin-flip line from neutral hydrogen. While those ex- 
periments are sensitive to the neutral gas distribution 
dominated by the IGM, CII will probe the onset of star 
formation and metal production in z ~ 6 to 8 galaxies. 

Recently it has also been proposed to use rotational 
lines of CO molecule to probe reionization (e.g., Gong 
et al. 2011; Carilli 2011; Lidz et al. 2011). CO studies 
have the advantage that redshift identification is facili- 
tated by multiple J transition lines, while with CII some 
confusion could result in the line identification with other 
atomic fine-structure and molecular lines at sub-mm and 
mm wavelengths. In comparison, and based on limited 
observational data, CII emission is expected to brighter 
in low mass galaxies, compared to the case of CO lumi- 
nosity. It is unclear if the galaxies present during reion- 
ization are analogous to low-redshift dwarf galaxies or 



the local ultra-luminous infrared galaxy population. At 
high redshifts, we expect most of the carbon to be in CII 
rather than CO, since the high-redshift galaxies may not 
host significant dust columns required to shield CO from 
dissociating UV photos. Given that a CO experiment 
to study reionization will involve an experiment at mid 
radio frequencies of 15 to 30 GHz, while CII will involve 
high frequencies, the two methods will be affected by dif- 
ferent systematics and foregrounds. Thus, a combined 
approach involving multiple probes of reionization, 21- 
cm, CO and CII, would provide the best avenue to study 
the z > 6 universe. 

In this paper we present an analytical calculation 
to predict the CII intensity intensity as a function of 
the redshift by considering spontaneo us, stimulated and 
collisional emiss ion processes ( jSuginohara et al.l 119991 : 
iBasu et aTll2QQ4f ). This intensity changes the brightness 
spectrum of the cosmic microwave background (CMB) at 
the frequency corresponding to the CII line (Basu et al. 
2004). In this paper, we focus on the CII flux from indi- 
vidual galaxies where the matter density is high and the 
collisional emission is the dominant process. As a check 
on our analytical calculations, we also consider results 
derived from numerical simulations to establish the CII 
intensity. The two approaches are generally consistent. 
We then consider the measurement of the CII intensity 
fluctuations resulting from the clustering of the galaxy 
distribution and sources that are present during and to- 
wards the end of reionization at z ~ 6 — 8. 

Experimentally, there are several challenges that one 
must overcome before the CII intensity fluctuations from 
high redshifts can be reliably established. First, higher J 
transitions of CO from dense molecular gas at lower red- 
shifts contaminate the CII line intensity measurements. 
In particular, one must account for all of CO (2-1) to 
CO (13- 12) emission lines from individual narrow redshift 
ranges in the foreground between and 6 when probing 
CII fluctuations at z > 6. To the extent to which a va- 
riety of existing; p redictio ns on the C O intens ity can b e 
trusted (|Gong et al.ll2QllULidz et aLll2011b ICariliai2011h , 
we generally find that the contamination is mostly below 
the level of the CII signal. Extending previous studies 



(IBasu et al. 2004: Visbal fc Lo eb 2010: IGonget al.l[20lTI : 



iLidz et al.ll2011[ ), we propose the cross-correlation of CII 
line intensity mapping and 21-cm fluctuations as a way to 
improve studies related to the epoch of reionization. To 
evaluate the detectability of the CII signal, we calculate 
the errors on the CII power spectra and CII-21cm cross 
correlation, respectively, based on the design of potential 
(sub-)milimeter surveys for the CII emission. For 21-cm, 
we consider the first generation experiment LOw Fre- 
quency ARra}{3 (LOFAR) survey as well as the improve- 
ment expected from a second generation experiment like 
the low-frequency extension to the Square Kilometer Ar- 
ray (SKA). 

The paper is organized as follows: in the next sec- 
tion, we derive the formulas to calculate the spin tem- 
perature of the CII line in the ISM of galaxies and the 
IGM. In Section 3, we calculate the mean CII intensity 
analytically and compare it with results derived from a 
simulation. We show the CII power spectrum in Sec- 
tion 4. In Section 5, we discuss the low-redshift contami- 

1 http://www.lofar.org/ 



nation from CO emission lines for CII intensity mapping 
at z > 6 and in Section 6 we propose a cross-correlation 
between CII and 21-cm line intensity measurements over 
the overlapping redshift ranges as a way to both distin- 
guish the CII signal from the CO emission and to improve 
overall constraints on reionization physics. In Section 7 
we outline the experimental parameters of a CII spec- 
tral mapping experiment designed to probe z ~ 6 to 8 
CII intensity fluctuations and discuss the detectability 
of the CII power spectrum and the CII and 21-cm cross- 
correlation when combined with LOFAR and SKA. We 
conclude with a summary of our results in Section 8. 
Throughout this paper, we assume the flat ACDM model 
with fi m = 0.27, Ct h = 0.046, a 8 = 0.81, n s = 0.96 and 
h = 0.71 (jKomatsuet al.ll2011l ). 

2. THE DERIVATION OF THE SPIN TEMPERATURE OF 
THE CII EMISSION 

The CII line intensity can be generated from carbon 
present in both the interstellar medium (ISM) of individ- 
ual galaxies and the diffuse intergalactic medium (IGM) 
in between galaxies. In the ISM of galaxies, CII is ex- 
pected to be in thermodynamic equilibrium with reso- 
nant scattering emission off of CMB photons and the 
collisional emission induced by gas in the galaxy. If the 
number density of gas particles n gas is greater than a crit- 
ical value n cr , the collisional excitation and de-excitation 
rate would exceed that of the radiative (spontaneous and 
stimulated) processes. 

Since the ionization potential of carbon is 11.26 eV, be- 
low 13.6 eV of hydrogen ionization, the surrounding gas 
of CII ions can be either neutral hydrogen or electrons. 
However, because the critical number density of electrons 
to trigger collisional excitation n™ is less than 100 cm -3 
while the critical number density of neutral hydrogen 
for collisional excitat ion ng is about 10 3 to 10 4 cm -3 
(jMalhotra et al.l l2001). the electrons can collide with CII 
ions mo re frequently than with HI, especially in ion- 
ized gas (jLehner et al. 2004; S uginohara et al.l fl999). For 
simplicity and not losing generality, we assume that the 
ISM is mostly ionized in individual galaxies. Then the 
CII line emission would main ly be excited by electrons 
in the ISM (| Suginohara et al.lll999l ). 

On the other hand, in the diffuse IGM in between 
galaxies CII line emission will be mainly due to ra- 
diative processes such as spontaneous emission, stim- 
ulated emission due to collisions with CMB photons, 
and a UV pumping effec t similar to the Ly-q coupling 
for the 21-cm emission ([W outhuvsenl 11952b iFieldl 119581 : 
iHernandez-Monteagudo et al.l [20061 ). We will focus on 
the CII emission from the ISM of galaxies first and then 
consider the signal from the IGM. The latter is found to 
be negligible. 

2.1. The CII spin temperature in the ISM of galaxies 

In the ISM of galaxies, the ratio of thermal equilibrium 
population of the upper level 2 Ps/2 and lower level 2 P\/2 
of CII fine structure line can be found by solving the 
statistical balance equation 



n u _ Biyly + n e C\ u 
n\ B U \I V + Ail + n e C u \ 

= — exp[-T* >u i/Ts >u i]. 
9\ 



(1) 



Here A u \ = 2.36 x 1 Q~ 6 s -1 is the spontan eous emission 
Einstein coefficient (jSuginohara et afl l999). B u \ and B\ u 
are the stimulated emission and absorption coefficients, 
respectively, I v = B[Tcmb(z)] is the intensity of CMB 
at z/ u i, n e is the number density of electrons, and C\ u 
and C u \ are the excitation and de-excitation collisional 
rates (in cm 3 s _1 ), respectively. Note that the UV pump- 
ing effect is neglected here and, as we discuss later, it 
should not affect the result unless the UV intensity inside 
the galaxy (/§y) is higher than 10 -15 erg/s/cm 2 /Hz/sr, 
which is about 1Q 6 times greater than the UV back- 
ground (jGiallongo et al.l 119971 : iHaiman et al.ll2QQQ[ ). We 
note that even if the radiative coupling from UV pump- 
ing made a small contribution to the spin temperature it 
would be in the same direction as the collisional coupling 
since the UV color temperature follows the gas temper- 
ature. 

The second line of EqJT] defines the excitation or spin 
temperature Ts, u \ of the CII line. The statistical weights 
are g u — 4 and g\ = 2, and T* ?u i = hv u \/kB — 91K is 
the equivalent temperature of the level transitio n. The 
excitation collisional rate C\ u can b e written as (|Spitzerl 
[19781 IOsterbrockl[T98l E[val 2008) 



Clu 



8.629 x 10~ 6 



7iu(T k e )exp 



-T*. 



lu 



(2) 



where T£ is kinetic temperature of the electron and 7i u is 
the effective collision strength, a dimensionless quantity. 
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Fig. 1. — The CII spin temperature T§ u and the T C mb in the 
ISM of galaxies as a function of the redsnift. To capture different 
ISM conditions within a large sample of galaxies we select several 
pairs of the electron kinetic temperature T^ and the number density 
n e for this calculation. We find that the CII spin temperature is 
almost constant with redshift, which indicates that the collisional 
process is dominant in the ISM of galaxies when compared to the 
radiative processes. 

To derive the spin temperature Ts, we can use Ein- 
stein relations g\B\ u = g u B u \ and A u \ = (2hv 3 /c 2 )B u \ 



to convert both stimulated emission coefficient B u \ and 
absorption coefficient B\ u in terms of the spontaneous 
emission coefficient A u \. Also, using the collisional bal- 
ance relation 



Clu _ 
Cul 

we can finally write 



^exp[-T, )Ul /T k e ], 



3i 



(3) 



T*. 



ul 



Ts,u] 



log< 



A ul (I uC y2hu s )^ 



n e C u iexp(-T* )U i/r^) J 

(4) 
Note that the de-excitation collisional rate C u \ is de- 
pendent on the Tfc , and using Eq. [2] and Eq. [3j we can 
calculate the de-excitation collisional rate C u \ for a fixed 
value of the electron kinetic temperature. Here we adopt 
the values of 7i u which are calculated by the R-matrix in 
Keenan et al. (1986), and they find 7i u = 1.58, 1.60 and 
2.11 at T£ = 10 2 , 10 3 and 10 4 K respectively. Putting all 
values together, we find C u i = 3.41 x 10" 7 , 1.09 x 10" 7 



and 4.55 x 10" 



for TS = 10 2 , 10 3 and 10 4 K, re- 



spectively. This is well consistent with the de-excitation 
rate o f 4.6 x 10 -8 cm 3 s _1 at 10 4 K given in iBasu et al.l 
(J200l . 

The deviation of the CII spin temperature Ts, u i rela- 
tive to Xcmb as a function of redshift for the electron 
dominated ISM of galaxies is shown in Fig. [TJ We choose 
several values of the electron kinetic temperature T^ and 
the number density n e to plot the spin temperature Ts 
as a function of the redsnift. We note that although 
the mean number density of electrons n™ ean can be very 
small (< 10 -2 cm -3 ) in a halo, the significant CII emis- 
sion in galaxies comes fro m dense gas clumps which have 
a much higher n e (Suginohar a et al.l fl999). Due to gas 
clumping the local n e can be much greater than n™ even 
though n™ ean < n™ . Thus we choose to assume several 
high n e values in our plots, but also show the case with 
two low values of n e = 1 and 10 cm -3 . These values are 
less than nf '. We find the Ts is almost constant and al- 
ways greater than Xcmb for < z < 10 in all these cases. 
In Eq. 21 it is easy to find that the Ts depends on the 
relative strength of the radiative (spontaneous and stim- 
ulated) and collisional processes. If the spontaneous and 
stimulated emission are dominant, we have Ts ~ Tcmb, 
while Ts ~ T£ if collisions are dominant. Given a fixed 
number density n e , the only variable that depends on 
redshift z in Eq. [4] is Iv{z), but the spin temperature is 
not strongly sensitive to it. This implies that the colli- 
sional process is dominant in the ISM of galaxies when 
compared to the resonant scattering off of CMB photons. 
As we discuss next, this result is not true for the emis - 
sion of the CII line in the diffuse IGM (|Basu et al.ll2004D . 
where Ts is much smaller and varies with redshift similar 
to Tcmb- 

2.2. The CII spin temperature in the diffuse IGM 

In the IGM, the collisional process becomes unimpor- 
tant since the number density of electrons and other el- 
ements are much s maller than in de nse regions within 
the ISM of galaxies (|Basu et al.ll2004J ). Also, the sponta- 
neous emission and the stimulated absorption and emis- 
sion by the CMB photons are considerable. We will also 
take into account UV pumping that can enhance CII 
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Fig. 2. — The CII spin temperature T]g n and T C mb in the diffuse 
IGM as a function of the redshift. We select several pairs of the 
UV color temperature Tuv an d the UV background intensity l}? v 
for this calculation. We find the spin temperature is close to Tcmb 
at high redshifts [z > 2). 

emission. This effect is si milar to the Wo uthuy sen- Field 
effect for the 21-cm line (| WouthuYipnl H952|; lFieldlll958b 
iHernandez-Monteagudo et al.l 120061 ) and, to the extent 
we are aware of, has not been discussed in the literature 
previously. 

At high redshifts the soft UV background at 1330 A 
generated by the first galaxies and quasars can pump the 
CII ions from the energy level 2s 2 2p 2 P\/2 to 2s2p 2 2 D 3 / 2 
at 1334.53A and 2s 2 2p 2 P 3/2 to 2s2p 2 2 D 3/2 at 1335.66A. 
Then this pumping effect can lead to the CII fine- 
structure transitions 2 D 3 / 2 ~^ 2 Pzii ~^ 2 Pi/i, which 



'3/2 — r ^3/2 — > ^1/2 > 

would mix the levels of the CII line at 157.74/im. The 
UV d e-excitation and excitation rates are given by ([Fieldl 

mm) 



and 



dUV 
ul 



pUV _ 

-flu 



9k 



9u En ^kn 



A 



ku 



9k 



A 



ku 



91 En ^kn ' 



c 2 4 UV 

2hv 3 



2hv 3 



(5) 



(6) 



where "k" stands 

18 c-l 



2.41 x 10* s- 1 and A ku 
coefficients of 2 Dq 



for the level 2 D 3/2 , #k = 4, A kY = 
4.76 x 10 T s _1 are the Einstein 
3/2 ^ 2 P1/2 and 2 ^3/2 ^ 2 P3/2 respec- 
tively (see the NIST website 0), and E n ^kn = A k i+A ku . 
Also, it is helpful to define the UV color temperature 
2~uv,ui m terms of the ratio of the UV de-excitation and 
excitation rates 



pUV 



9u r 

— exp- 

9\ 



■^•,ul/^UV,ul]- 



(7) 



Then, following a similar derivation as to the case of 

2 http://physics.nist.gov/PhysRefData/ASD/lines_form.html 



spin temperature in the ISM of galaxies, we find that the 
CII spin temperature in the diffuse IGM is 



Ts,u 



log< 



A ul [l + (I„c 2 /2hv s )„J 



T 



A ul (I u c 2 /2h^) Uul + Puv e xp(-T^ ul /T uv ) J 

(8) 
In Fig. [2j we show the departure of Ts from Tcmb as 
a function of the redshift. We see that the T<p M does 
not change if Tuv > 10 3 K when fixing I^ v , and it is 
close to Tcmb at high redshifts (z > 2). Together with 
Eq. IHlthis implies that A u \ ^> P^j v , i.e. the spontaneous 
and stimulated emission are much greater than the UV 
pumping effect at high redshifts. By comparing Fig. [2] 
with Fig. [TJ we can see that the CII spin temperature 
in the ISM of galaxies is much larger than that of the 
diffuse IGM, and T<p M is quite close to T CM b at high 
redshift while T| a is larger than Tcmb- 

3. THE CALCULATION OF THE CII MEAN INTENSITY 

To estbalish the overall intensity of the CII line emis- 
sion we calculate the distortion of the CMB spectrum 
AI U and take into account the mean intensity of the CII 
line emitted by the ISM galaxies and the IGM. Consider- 
ing the spontaneous, stimulated emission and absorption, 
and the expansion of the Universe, we write the radiative 
transfer equation as 



dl u (z) 
ds 



■ j y {z) - a v (z)I v 



3^, 



(9) 



where ds is the line element along the line of sight, and 
H{z) is the Hubble parameter. The spontaneous emis- 
sion and absorption coefficients are 



3v{z) 



hv u \ 

47T 



n u (z)A u i<j)(v) 



and 



OL y {z) 



hv n \ 

47T 



(j){v)(n\B\ u -n u B u i) 



respectively, where 4>{v) is the line profile function which 
can be set as delta function 4>{v') = 8{y' — v) if the 
thermal broadening o r velocity is much smal l er than the 
frequency resolution (Sugi nohara et aflll999( : iBasu et all 
2004). Integrating Eq|9] along the line of sight then gives 






j v (z) - a v (z)I v 



(1 
3v{z) - 



Vzf 

OL v {z)I v 



H(z)(l + z) 4 



ds 



cdz. 



(10) 



Using the relation of the line profile and the redshift 
4>(v') = 0[^o(l + z')\ = [(1 + z)/v\ 5{z' — z) we can obtain 
the integrated intensity of the CII emission lines at z as 



AI U 



tonm + iy ** f ™ n ™ {z) x (ID 



^exp(-T*, u i/r SiUl ) 

f grd 



exp(T*, ul /T s , ul ) 



(2hi^/c 2 I v ) Vu 



where f^ u is the fraction of CII ions at the ground level 
2 P\j2- If Ts ^> T*(> Tcmb), which is the usual case in 



the galaxies, then ex.p(±T* yU \/Ts yU i) ~ 1, and we can find 



AT 



he 



9u 



47rH(z)(l + z) 3 gi 



A ui /cii n cn(^). 



(12) 



If Ts is much larger than T* jU i (e.g. the cases in galax- 
ies), /qjj ~ 1/3 is a good approximation, and when 

Ts <C T* jU i we can set /§jj c± 1 (e.g. the case of the 
IGM). The ncn(^) is the number density of the CII ions 
at z, which has to be carefully estimated for both the 
galaxy and IGM cases. 

To calculate the CII intensity, we now need an estimate 
of the number density of the CII ions. For the IGM case, 
it is easy to calculate if we know the metallicity evolu- 
tion and the average baryonic density as a function of z. 
However, for the ISM of galaxies, we also have to find the 
fraction of the CII ions which exceed the critical density 
to trigger the dominant collisional emission. In the next 
section, we will first evaluate the ncn theoretically for 
both the ISM and IGM cases, and then calculate ncn 
using data from a simulation and compare it to the an- 
alytical result as a check on the consistency. We also 
provide an order of magnitude estimate on the z ~ 6 
to 8 CII mean intensity using scaling relations such as 
those involving star-formation rate and the total far-IR 
luminosity of galaxies. 

3.1. The analytic estimation 
We start by writing 

ncn(z) = fcn(z)Z c n gas (z), (13) 

where Z c = X C Z% , Z% = 3.7 x 1(T 4 is the solar 
carbon abundance (jBasu et all I2QQ4I ) , and we assume 
X c = 1 in th e galaxy and X c = 10~ 2 in the IGM 
dSavagliol 119971 : lAguirre fc Schavd 120051 Kramer et all 
2010). The fcu(z) = Z/Z Q (z) is the global metallicity at 
different redshift, and we use an approximated relation 
Z/Zq(z) = 10 -0,42 in our calculation, which assumes the 
present mean metallicity of the Universe is equal to the 
solar metallicity (note that this assumption may over- 
estimate the metallicity at low redshifts) and the car- 
bon atoms are totally ionized to be CII. This relation is 
consistent with the observational data from the damped 
Ly-a absorbers (DLAs) metallicity measureme nts, which 
cover s the redshift range from 0.09 to 3.9 (jKulkarni et al.l 
2005 ), and also matches a previous theoretical estim ation 
(jPei et al.lll995l fl99l iMalanev fc Chaboverl 1 1996ft . 

The n gas (z) in Eq. [13] is the mean number den- 
sity of the gas, which, for the IGM case, is just the 
baryon density, fih(z) = 2 x 10 _T (1 + z) 3 cm -3 . For 
the ISM of galaxies, we use n gas (z) = fg&sfgl s nb(z), 
where f™ is the fraction of the gas that are present 
in d ense environments of th e ISM and satisfies n gas > 
n cr (jSuginohara et all Il999f ). The value of f g l s de- 
pends on the gas clumpin g within the ISM of galax- 
ies and the Jeans ma ss (jFukugita fc Kawasakil Il994t 
iSuginohara et al.lll999l ). We find f™ s is about 0.1 at 
z = 6 and decreases slowly a t higher redshifts in our 
simulation (|Santos et al.l l2Q10h . For simplicity, we just 
take /" = 0.1 to be the case for all galaxies indepen- 
dent of the redshift. As is clear, this parameter is the 
least uncertain of the calculation we present here. Ob- 



servations with ALMA and other facilities of Herschel 
galaxy samples will eventually narrow down this value. 

The fraction of the "hot" gas (T > 10 3 K) in halos 
/g°s is also included here, since the main contribution 

of CII emission comes from the gas with T > 10 3 K 
(see Fig. [3|). We find the fraction is around 0.3 and re- 
main s constant from z = 8 to 6 in the De Lucia cata- 
log (jDe Lucia fc B laizot 2007), which is used to derive 
the CII number density from the simulation in the next 
subsection. As we are computing the mean intensity ex- 
pected in a cosmological survey, n gas is the number den- 
sity for galaxies in a large space volume instead of one 
individual galaxy and thus we can still use n^ for the 
ISM of galaxies. 
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Fig. 3. — The mean intensity of CII emission line from the ISM 
gas in galaxies as a function of the redshift. We find that the 
electron collisional emission saturates when n e > 10 2 cm -3 and 
T^ > 10 3 K. 

In Fig. [3j we plot an analytical estimate of intensity 
of CII emission line from the ISM of galaxies as a func- 
tion of redshift. We select the same pairs of T£ and n e 
as in the Fig. [TJ We see that the intensity is practically 
independent of the electron number density and the tem- 
perature when n e > 10 2 cm -3 and T£ > 10 3 K, e.g., the 
signal is only seen in emission and is essentially propor- 
tional to the CII density (see Eq. [T2]) . Even for n e = 1 
and 10 cm -3 which is less than the n^ 1 *, we still find a 
significant intensity for the CII emission from the ISM of 
galaxies. 

In Fig. [H we show the same for carbon in the diffuse 
IGM. We find the I™ from the IGM is much smaller 

than that from galaxies, and ^cn/^cn M ~ 10 4 f° r a ^ 
cases we consider at all redshifts. This is because the 
CII spin temperature and the CII abundance in the ISM 
of galaxies are much larger than that in the diffuse IGM. 
Thus the CII emission from IGM can be safely neglected, 
and hereafter we will just take into account the CII emis- 
sion from the ISM of galaxies when discussing intensity 



6 



fluctuations. 

Note the line intensity measurements of individual 
galaxies are generally described with PDR models using 
the number density and the UV intensity within the ISM, 
inste ad of number density and temperature as we use 
here dMalhotra et all 12001 ). We depart from the PDR 
modeling approach as we are considering global models 
that are appropriate for cosmological observations and 
are not attempting to model the properties of line emis- 
sion from individual galaxies. It is likely that future work 
will need to improve our model by making use of a more 
accurate description of the stellar mass function and the 
star-formation history of each of the galaxies in a large 
cosmological simulation by performing calculations to de- 
termine T e given the UV intensity field dominated by 
the massive, bright stars. When making predictions re- 
lated to the CII intensity power spectrum (Section 4), 
we take as default values T e = 1000K and n e = 100 
c m" 3 . These values are f ully consistent with the findings 
of iMalhotra et all (j2QQlh . where they find for 60 normal 
star-forming galaxies, the values of T ~ 225 to 1400K 
and n ~ 10 2 to 10 4 cm -3 . 

3.2. Intensity estimated from numerical simulations 

Since the evolution of the metallicity and the critical 
fraction of the gas f™ s are hard to estimate, the num- 
ber density of the CII ions ncn is not well determined 
analytically for the case involving CII emission from the 
ISM of galaxies. So we use results derived from a numer- 
ical simulation to check the analytical estimation. To do 
this we first derive the expected CII mass, Men, in a 
halo with mass M at a given redshift. Then, by inte- 
grating over all possible halo masses in a given volume, 
we estimate the CII mass for that same volume. 

To find Mcn(z) with simulations, we make use of the 
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Fig. 4. — The mean intensity of CII emission line from carbon in 
the IGM as a function of the redshift. The intensity here is much 
smaller than that of ISM in galaxies since the spin temperature is 
close to the CMB temperature and the gas density is much lower 
than galaxy ISM. 
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Fig. 5. — The mass in metals Mz as a function of the halo 
mass M at z = 8. The solid lines shows the mean relation (blue 
center line) and ±1<7 relation (green lines). The red dots denote 
the mean value of the scattering (gray region) when binned to 150 
logarithmic intervals in halo mass. 



galaxy catalog from iDe Lucia fc Blaizotl ([20071 ). This 
catalog is obtained by post-processing the Millennium 
dark matter simul ation with sem i -ana lytical models of 
galaxy formation ([Springel et al.l l2005f h and has the 
metal content in each of four different components: stars, 
cold gas, hot gas and an ejected component. In this cal- 
culation we will assume that the CII emission comes from 
the hot gas component which has a te mperature of about 
IO 5 - 1 7 K. However, according to Obresch kow et all 
(2009a), some of the gas in galactic disks which is con- 
sidered as cold gas (T « IO 2 — IO 3 K) in the De Lucia 
catalog, should actually be warm (T w 10 4 K) and ion- 
ized. So we will al so consider the case where 10% of this 
cold gas from the IDe Lucia fc Blaizotl ([20071 ) simulation 
is recategorized as warm and thus contribute to the CII 
luminosity. 

In Fig. [5] we plot the mass in metals Mz as a function 
of the halo mass M for the hot gas only case. There is 
some scatter, specially at the low mass end, but note that 
if the experimental volume is large enough, we will be ob- 
serving a large number density of galaxies such that the 
total Mz should approach the average. The average rela- 
tion between Mz and halo mass M can be parameterized 
in the form M Z (M) = M (M/M C ) 6 (1 + M/M c )~ d . At 
z = 6, z = 7 andz = 8, these parameters take the values 
M = 1.1 x IO 9 , 1.0 x IO 9 and 1.6 x IO 9 , M c = 3.5 x 10 11 , 
3.5 x 10 11 and 3.7 x 10 11 , b = 3.6, 3.4 and 3.4, and 
d = 3.25, 3.1 and 3.6, respectively. 

At high ISM temperatures with T > 100K we assume 
carbon is ionized, so we have Mcii=/q Mz where f^ = 
0.21 is the carbon fraction of the mass in metals in the 
Sun (Arnett 1996). By taking into account the expected 
number of halos, we can obtain ncn as 



»ssc 



0=/ 



dM dn_M cu (M,z)^ 



dM 



m. 



where m c is the atomic carbon mass, M is the halo mass 
and dn/dM is the halo mass function (Sheth & Tormen 



1999). The integration is made from M min = 10 8 M /h 
to M max = 10 13 M /h, and M m i n is the minimum mass 
of the dark matter halo that can host galaxies. Note that 
the result is insensitive to the exact value of M max since 
halos with M > 10 13 are rare due to the exponential term 
in the halo mass function. From the previous section we 
can safely assume that the spin temperature is saturated 
(e.g. the CII is only seen in emission) so that we now 
have all the ingredients to calculate the signal using the 
simulation values. 




Fig. 6. — The CII proper number density as a function of redshift. 
The red solid line is derived from our simulation with just hot gas 
contributing to Lqii, the yellow dotted line is our analytic result, 
and the blue spots are obtained from the simulation using both the 
hot gas and the warm gas (10% of the cold gas in the galaxies; see 
text for details). 
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Fig. 7. — The mean intensity of CII emiss i on lin e from the out- 
puts obtained from the De Lucia & Blaizot (2007) simulation at 
several different redshifts. The yellow dotted line is our analytic 
result, and the blue spots are obtained from the simulation with 
the hot and warm gas contributing to the Lqii- The yellow line and 
the blue spot were calculated with n e = 10 3 cm -3 and T^ = 10 4 K. 
The other lines are derived from the simulation assuming only the 
hot gas contributes to the Leu- 



In Fig. [6j we show the number density of the CII ions 
n cu( z ) usm g the simulation result. We find our analytic 
result (black dashed line) is well consistent with that de- 
rived from the simulation especially at high redshift. The 
departure at low redshift in our analytical calculation is 
reasonable since the the hot gas fraction f^ should be 
higher at low redshifts because of the higher star forma- 
tion rate. 

Assuming only the hot gas component of 
iDe Lucia fc Blaizotl (J2QQ7D contribute to the CII 
line luminosity, Leu, our two calculational methods 
differ by a factor of a few; however, after adding the 
warm component, the analytical result of the previous 
section match the numerical values at z > 6 where we 
are especially interested. The difference is again same 
for the mean intensity of the CII emission shown in 
Fig. [7J Again, we find that the CII emission is saturates 
when n e > 10 2 cm -3 and T£ > 10 3 K. For higher values, 
CII intensity is almost independent of n e and T£. 

3.3. Intensity estimated from observed scaling relations 

At z ~ 7, our calcu lations involving either analytical or 
through lDe Lucia fc B laizot (2007) simulations suggest a 
mean intensity of about 100 to 500 Jy/sr with a preferred 
value around 200 Jy/sr. We can use an approach based 
on obse rved scaling relations, similar to the approach 
used by Carilli (2011) to estimate the mean CO (1-0) in- 
tensity during reionization, to estimate the mean CII 
intensity. We begin with the star-formation rate from 
z ~ 6 to 8. While estimates exist in the literature from 
Lyma n break galaxy (LBG) dropouts (|Bouwens et al.l 
2008), such estimates only allow a lower estimate of the 
star-formation rate as the luminosity functions are lim- 
ited to the bright-end of galaxies and the slope at the 
low-end is largely unknown and could be steeper than 
indicated by existing measurements. An independent es- 
timate of the star-formation rate duri ng reionization in - 
volve the use of gamma-ray bursts (iKistler et al. 2009). 
Together, LBG and GBR star-formation rates lead to a 
range of 0.01 to 0.1 M yr _1 Mpc~ 3 at z ~ 7. 

We can convert this SFR to the CII luminosity of 
~ 3 x 10 40 to ~ 3 x 10 41 ergs s" 1 usi ng the observed 
scaling relation of IDe Looze et all (|2011[ ) for low-redshift 
galaxies when averaging over 100 Mpc 3 volumes. We also 
get a result consistent with this estimate when convert- 
ing the SFR to total integrated FIR luminosity and then 
assuming that 10~ 2 to 10 -3 of the Lfir appears in CII 
iStacev et al.l (|201Q[ ) , again consistent with z ~ 2 redshift 
galaxy samples. Once again we note that our estimate is 
uncertain if either the SFR to CII luminosity calibration 
evolves with redshift or the CII to FIR luminosity ratio 
evolves with redshift. The above range with an order of 
magnitude uncertainty in the SFR, and subject to uncer- 
tain evolution in observed scaling relations from z ~ 2 
to 7, corresponds to an intensity of 40 to 400 Jy/sr at 
z ~ 7. This range is consistent with our independent 
estimate, but could be improved in the near future with 
CII and continuum measurements of high redshift galaxy 
samples with ALMA. 

4. THE CII INTENSITY POWER SPECTRUM 

The CII intensity we have calculated is just the mean 
intensity, so in this Section we will discuss spatial varia- 
tions in the intensity and consider the power spectrum of 



the CII emission line as a measure of it. This power spec- 
trum captures the underlying matter distribution and if 
CII line intensity fluctuations can be mapped at z > 6, 
the line intensity power spectrum can be used to probe 
the spatial distribution of galaxies present during reion- 
ization. 

Since the CII emission from the ISM of galaxies will 
naturally trace the underlying cosmic matter density 
field, we can write the CII line intensity fluctuations due 
to galaxy clustering as 



CII; 



sir = bcuit; n 6( x ) 



(15) 



where I^ u is the mean intensity of the CII emission line 
from the last section, 6(x) is the matter over-density 
at the location x, and bcu is the average galaxy bi as 
weighted by CII luminosity (see e.g. Gon g et al.ll201l[ ). 
Following Eq.[l4]and taking into account that the fluc- 
tuations in the halo number density will be a biased 
tra cer of the dark matter , the average bias can be written 
as (jVisbal fc Loebll20T0h 



bcn(z) = 



fM:::d M 7M M cu(M)b(z,M) 

: j ( 16 ) 

fdn 71//-/ fl/A 
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aiV1 dM 



diLMcn(M) 



where 6(z, M) is the halo bias and d n/dM is the halo 
mass function ([Sheth fc Tormen| [l999). We take M min = 
10 8 M /h and M max = 10 13 M /h. Then we can obtain 
the clustering power spectrum of the CII emission line 



P^i(z,k)=b z CII IcnPss(z,k), 



(17) 



where Pss(z, k) is the matter power spectrum. 

The shot-noise power spectrum, due to discretization 
of the galaxies, is also considered here. It can be written 
as ([Gong et al.ll2Qllh 



Dshot 
CII 
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J M mi: 



^m{^ )Dl 



(18) 



where D\, is the luminosity distance, .Da is the comov- 
ing angular diameter distance and y(z) = dx/dv = 
Acn(l + z) 2 1 'H(z), where \ is the comoving distance, v is 
the observed frequency, Acn is the rest frame wavelength 
of the CII line. The Leu is the CII luminosity which can 
be derived from the i^? n , and we assume Lcii(M, z) = 
B x M z {M,z) [£©] and finally find B = 100.63, 100.36, 
and 100.17 at z = 6, 7, and 8 respectively. The total CII 
power spectrum is P&\(z, k) = Pgfi 8 (z, k ) + p c\\\ z )- 

In Fig. El as an example we plot the clustering, shot- 
noise and total power spectrum of the CII emission at 
z = 6, z = 7 an d z — 8. Using the Halofit code 
(Smi th et al.l 120031 ). we calculate the non-linear matter 
power spectrum. In order to calculate the intensity of 



the CII line we use T£ = 



10 J K and n e = 10 2 cnT 



which are a possible representat ion of the conditions at 
which CII emits in the galaxies (Mal hotra et al.l [200 If ). 

For comparison, the CII power spectrum estimated 
by other works are also shown here. The blue dash- 
dotted line denotes the P^\ derived from the rela- 
tion Lcn/Lco(i-o) — 10 4 ( IBreuck et al.l [20111). and 
the Lqo(i-q) is from the calculation of IVisbal fc Loebl 
(|2010f ). The magenta long dashed line is the P^\ evalu- 



ated v ia the CII luminosity derived from lVisbal fc Loebl 
(J2010l ). 

As expected the CII power spectrum is larger than 
the C O(l-O) power spectrum calculated in iGong et all 
(201 l]). N o te th at the CO intensities predicted by 
IGong et al.l (J2011I ) should be corrected by a factor of 
1/(1 + z) which comes from a missing conversion factor 
between the CO flux obtained from Obreschkow et al. 
(|2009b[ ) and the actual CO(1-0) luminosity. Although 
this correction further reduces the CO signal, we point 
out that the model we assumed for the CO luminosity 
as a function of the halo mass generates a steep depen- 
dence at the low mass end. If we use a less steep model at 
the low end of halo masses (M < 10 10 M /h), where the 
simulation has a large scatter, such as ^co(i-o) <* ^haio, 
then the CO signal can increase by a factor of a few, par- 
tially compensating the 1/(1 + z) correction. 

Comparing the CO (1-0) power spectrum (using the 
1/(1 + z) correction but keeping the luminosity model 
used in IGong et al.l (J2011f )) with the CII one with its 
density estimated using only the hot gas component 
from the simulation, implies a luminosity relation of 
£cn/£co(i-o) ~4x 10 4 when z ~ 6 reducing to below 
10 4 when z ~ 8. The observational value for this rela- 
tion as valid for < z < 2 is ^ch/^co(i- o) ~ 4100 as 
obtained by iStacev et al.l (jl991f ) (see also iStacev et al.l 
2010). However detections of CII emission from star- 
forming galaxies at z > 6 is almost non-existent and 
there is a strong possibility for an evolution with red- 
shift for this relation. For the handful of galaxies stud- 
ied, the possibility for evolution is supported by the 
£cn/£co(i-o) - 10 4 value derived for LESS J033229.4- 
275619 at z = 4.76 (jBreuck et al.ll2011f ). which is one of 
the highest redshift CII detections in a sub-mm selected, 
star- formation dominated galaxy by far, though higher 
redshift detections are soon expected as part of Herschel 
followup campaigns. In Fig. [8] we show that this observa- 
tional ratio (blue dash-dotted line) is consistent with our 
direct CII prediction from the simulation with the hot 
gas case (red solid line) around la level (green dashed 
lines, calculated using the la error of the Mz-M rela- 
tion). The noise power spectrum and the error bars are 
estimated for an assumed (sub-)milimeter survey, which 
we will discuss in some deta il in Section 7. 

The luminosity formula in lVisbal fc Loebl (J201Q[ ) based 
on the observed sub-mm and mm emission line ratios 
(see below for CO version) can also be used to calcu- 
late the the CII power spectrum as a function of the 
redshift, but it leads to a result (long dashed magenta 
line) that is smaller than that estimated by the Lcn- 
Lco(i-o) relation (blue dash-dotted line) and our CII 
power spectrum. This is effectively due to a difference 
in the calibration with CO (1-0) luminosity from M82 
and CII luminosity from a sample of low-redshift galaxies 
(Mal hotra et al.l 12001). The ratio ^cn/^co(i-o) can be 
as low a s 1.6 xlO 3 for l umin ous low-redshift galaxies such 
ULIRG aStacev et al.l (j2010f ). This sets a lower bound on 
the Leu estimation if we assume that z > 6 galaxies are 
analogous to local ULIRGs. 

In Fig. 9, we plot the inte nsity maps of the CI I emis- 
sion at z = 6, 7 and 8 using iDe Lucia fc Blaizotl ([2007) 
simulations to calculate the CII line intensities to show 
what the sky looks like at each of these redshifts when 
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Fig. 8. — The clustering, shot-noise and total power spectrum of CII emission line at z = 6, z = 7 and z = 8. The red solid, dashed and 
dotted lines denote the CII total, clustering and shot-noise power spectrum respectively. The Pq{{ s is estimated from the simulation with 
only the hot gas contributing to Leu, and we assuming that T^ = 10 3 K and n e = 10 2 cm -3 here. The green dashed lines are the la error 
of the CII power spectrum which are derived from the la error of M z — M relation in Fig. [5] The error bars and noise power spectrum (red 
long-dashed line) in the top left and top right panels are est imated with lm and 10m aperture for CII line respectively. The magenta long 
dashed line is derived from the CII luminosity estimated by Visbal & Loeb (2010). The blue dash-dotted line is estimated via the relation 
^Cll/^CO(1-0) — 10 4 - The bottom left and bottom right panels are the same plots but with the error bars on the total power spectrum. 



mapped with CII. The emission can be traced to individ- 
ual galaxies. Each of these maps span 3 degrees in each 
direction and the color bar to the right show the inten- 
sity scaling in units of Jy/sr. Note the large reduction in 
the maximum intensity from z = 6 to z = 8 due to the 
decrease in the overall abundance of metal content in the 
ISM of galaxies at high redshifts. 

5. CO CONTAMINATION TO CII LINE INTENSITY 
VARIATIONS 

When attempting to observe the z > 6 CII line in- 
tensity variations, the same observations will become 
sensitive to variations in the intensity of other emission 



lines along the line of sight. In particular low-redshift 
variations will be imprinted by a variety of CO transi- 
tions. This foreground line contamination distorts the 
CII signal and introduces additional corrections to the 
measured power spectrum, which will not be due to CII 
alone. 

As an example here we will focus on the contamination 
from the CO lines as we expect those intensities to be of 
the same order of magnitude as the CII line intensity. 
Given the difference in the CII rest frequency and those 
of the CO spectrum for each of the transitions, for a given 
redshift targeting CII line emission, the frequency of ob- 
servations and the bandwidth of the observations will 
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Fig. 9. — The sky maps of the CII emissi on at z = 6, 7 and 8 
from top to bottom, respectively using De Lucia & Blaizot (2007) 
simulations to calculate the CII line intensities. Each of these maps 
span 3 degrees in each direction and the color bar to the right show 
the intensity scaling in units of Jy/sr. 
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Fig. 10. — The contaminated CII total power spectrum P^^ 
(black solid lines) at z = 6, z = 7 and z = 8. The CII total power 
spectrum PqJj is calculated from the hot gas in the simulation (red 
solid line) (note that the PqI^ and Pq{\ are almost overlapped 
at z = 6 and z = 7 because of the relatively smaller CO line 
contamination). The other long dashed and green dotted line s 
are calculated with the ^co(m-n) given in lVisbal fc L oeb (2010). 
The green dotted lines from upper to lower are CO (7-6), CO (8-7) 
... CO(13-12) respectively. 

correspond to a certain low redshift range in which con- 
taminating CO contribution will be present . Here we cal- 
culate the CO power spectrum following Visbal & Loeb 
(2010) with a CO luminosit y formula as the fun ction of 
the redshift and halo mass ([Visbal fc Loebl 2010) 
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-CO(m-n) v..„, .,, \ 3 g x 1Q 6 y ^ 1O 1O M0 



M 



(19) 



"£©? 



l + z\ 3/2 /* 

7 / Cduty 

where Rn ne is the ratio of star formation rate and the 
luminosity for a emission line, /* = 0.1 is the fraction 
of gas in a halo that can form stars, and eduty = 0.1 is 
the duty cy cle which is canceled w hen computing the CO 
intensity in iVisbal fc Loebl (|2010l ) . 

We find this formula has some deviations from th e re- 
sults of our previous simulations (jGong et al.l l2011). but 
it is still a good approximation when considering the halo 
mass range 10 8 < M < 10 13 M /h in which we perform 
the calculation. The advantage of this formula is that 
we can calculate the luminosity of the CO lines at an 
arbitrary halo mass and redshift. Also, note that take 
assume Leo ~ M instead of Leo ~ Mqo when calculat- 
ing the large-scale structure bias factor of CO emitting 
galaxies. 

We show how the CO contamination for CII inten- 
sity measurements from a variety of redshift ranges for 
J = 2 — 1 to J = 13 — 12 transitions of CO combine 
in Fig. 10. The frequency of the CII line which is emit- 
ted at high redshift, e.g. z = 7, will be redshifted to 
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Fig. 11. — The power spectrum of 21-cm emission at z = 6, 7 and 8. The error bars in the left panel are estimated with a setup similar 
to LOFAR while the ones in the right panel are from a setup similar to SKA. 

mapping experiment be considered with another tracer 
of the same high redshift universe. In particular, low- 
frequency radio interferometers now target the z > 6 
universe by probing the neutral hydrogen distribution 
via the 21-cm spin- flip transition. Thus, we consider 
the cross correlation of the CII line and the 21-cm emis- 
sion at the same redshift to eliminate the low-redshift 
contamination, foregrounds, and to obtain a combined 
probe of the high redshift universe with two techniques 
that complement each other. We expect a strong cor- 
relation between the CII and 21-cm lines because they 
both trace the same underlying density field and such 
a cross-correlation will be insensitive to CO lines since 
they are from different redshift which depart far away 
from each other. There could still be lower order effects, 
such as due to radio foregrounds. For example, the same 
galaxies that are bright in CO could also harbour AGNs 
that are bright in radio and be present as residual point 
sources in the 21 cm data. Then the cross-correlation 
will contain the joint signal related to low redshift CO 
emission and the residual point source radio flux from 
21-cm observations. The issue of foregrounds and fore- 
ground removal for both CII alone and joint CII and 21- 
cm studies are beyond the scope of this paper. We plan 
to return to such a topic in an upcoming paper (Silva et 
al. in preparation). 

Here we calculate the power spectra for the CII- 
21cm correlation using the correlation between the mat- 
ter density field and the 21cm brightness temperature 
obtai ned from a simula tion made using the Simfast21 
code (Sa ntos et al.11201 0). with the further modifications 
described in (S antos et al.l l2011h to take into account 
the unresolved halos. This code uses a semi- numerical 
scheme in order to simulate the 21cm signal from the 
Reionization Epoch. The simulation generated has boxes 
with a resolution of 1800 3 cells and a size of L=lGpc. 
With an ionizing efficiency of 6.5 we have obtained the 



^o = ^cn/(l + z) — 237 GHz at present. As shown 
in Fig. [10] the main contamination is from the first five 
CO lines, CO(2-l) to CO(6-5), and the contamination 
of the CO lines provide about 2% and 30% to the to- 
tal intensity power spectrum at z = 7 and z = 8, re- 
spectively, for large scales (k < 1 h/Mpc). The black 
solid lines shown in the plot are the contaminated CII 
total power spectrum for the hot gas case of our simula- 
tion, which is the sum of the CII total power spectrum 
^cn(^) ( re d solid line) and the CO total power spec- 
trum P^Q^k) (dashed and dotted lines). At small scales 
(k > 1 h/Mpc), the shot noise becomes the dominant 
component in the power spectrum, and we find the shot 
noise of the CII emission is generally greater than that of 
the CO emission. We note that this result is likely sub- 
ject to assumptions on the CO luminosity caused by the 
CO luminosity we use, since the CO luminosity and es- 
pecially the CO luminosity-halo mass relati on. We find a 
weaker dependence on the halo mass than lVisbal fc Loebl 
(2010) in our CII emission model. 

The other atomic fine- structure emission lines at longer 
wavelengths than CII, for example Nil at 205 fim and 
CI with 370 and 609 fim can also contaminate CII emis- 
sion. The ratios of Lqu/L^u and Lcii /Lgi are about 
8 ~ 10 according to th e observations ([Malhotra et al.l 
120011 ; lOberst et aLll2006f h while their luminosities are ei- 
ther comparable or slightly lower than the CO luminos- 
ity. A CII study of reionization will also involve cross- 
correlation studies at low redshifts to untangle various 
line contaminations. In a future paper we will return to 
such a study using a numerical setup of a CII experiment 
of the type we propose here. 

6. CROSS-CORRELATION STUDIES BETWEEN CII AND 
21-CM OBSERVATIONS 

Since the above described CO contamination lines 
come from different redshifts it is necessary that any CII 



12 



io 5 



io 3 



10 2 



w IO 1 =* 



10° 



10- 



io- ; 



i i i m i 1 — i i i i m i 1 — I - n 



z = 6 T 



i* H *hi>. 



lm with LOFAR 



i i i i m i 1 — i i i i m i 1 — |— r 



z = 7 




rmrq | 

z = 8 



10" 1 10° 

k (h/Mpc) 



IO" 1 10° 

k (h/Mpc) 



IO" 1 10° 

k (h/Mpc) 




IO" 1 10° 

k (h/Mpc) 



O" 1 10° 

k (h/Mpc) 



Fig. 12. — The cross power spectrum of the CII and the 21-cm emission line at z = 6, z = 7 and z = 8. The red dashed lines denote 
the negative correlations while the red solid lines denote positive correlations. The la errors are also shown in thin red lines. Here we just 
consider the 21-cm signal from the IGM, since the 21-cm signal from neutral gas in galaxies is about 10 -4 of the IGM signal (Gong et al. 
2011). The error bars in the left panel are estimated using the (sub-) milimeter survey with aim aperture for CII line with a setup similar 
to LOFAR for 21-cm measurements, while 10 m aperture for CII line with a setup similar to SKA for 21-cm measurements in the right 
panel. 

mean neutral fraction fraction of 0.05, 0.35, and 0.62 and 
an average brightness temperature of 0.63 mK, 6.44 mK, 
and 14.41 mK for z = 6, z = 7, and z = 8 respectively. 
The power spectrum of 21-cm emission is also shown in 
Fig. EJ the blue dashed error bars are estimated from 
the LOFAR and the red solid ones are from the SKA (see 
Table 2 for experimental parameters). 

In Fig. [I2j we show the cross power spectrum of the 
CII and 21-cm emission line (red thick lines) and la un- 
certainty (red thin lines) at z = 6, z = 7 and z = 8. The 
error bars of the cross power spectrum are obtained by 
the assumed milimeter spectrometeric survey with 1 m 
and 10m aperture for CII line and LOFAR (left panel) 
and SKA (right panel) for 21-cm emission. Note that 
the correlation is negative on large scales (small k) when 
the correlation between the ionization fraction and the 
matter density dominates and positive on small scales 
when the matter density auto-correlation is dominating. 
This can be seen by looking at the expression for the 
cross-power spectrum, which to linear order is given by 



^cii,hi(2, k) oc 4/3(1 - Xi)P 66 - XiP 6x6 , 



(20) 



where x- x is the average ionization fraction and Ps x s is the 
cross power spectrum of the ionized fraction and the dark 
matter. To show the cross-correlation illustratively, we 
also plot the cross-correlation coefficient at z = 6, z = 7 
and z = 8 in Fig. [13j which is estimated by rcn,Hi(fc) = 

Pcn,m(k)/ \/Pcu(k)Pm(k) . The typical ionized bubble 
size can be reflected by the scale that rcn,Hi(fc) rises up 
from —1 towards zero. As indicated by the r*cn,Hi(fc) at 
z = 6, 7 and 8 in the top panels of Fig. [13j the ionized 
bubble size at z = 6 is greater than that at z = 7 and 8 
which is consistent with the current reionization scenario. 



7. OUTLINE OF A CII INTENSITY MAPPING 
EXPERIMENT 

We now discuss the requirements on an instrument de- 
signed to measure the CII line intensity variations and 
the power spectrum at z > 6. Here we provide a sketch 
of a possible instrument design including calculated noise 
requirements based on current technology; a detailed de- 
sign study will be left to future work as our concept is 
further developed. 

An experiment designed to statistically measure the 
CII transition at high redshift requires a combination of 
a large survey area and low spectral resolution; because 
of this, the use of an interfe r omete r array such as the 
one suggested in iGong et al.l (J2011f ) for CO seems im- 
practical. Instead we opt for a single aperture scanning 
experiment employing a diffraction grating and bolome- 
ter array in order to provide large throughput. A set of 
useful equations is outlined in the Appendix to convert 
from bolometer noise to the noise of the actual intensity 
power spectrum measurement. 

For proper sampling of the CII signal on cosmologi- 
cal scales and cross-correlation with 21-cm experiments, 
we would require at minimum a survey area of 16 deg 2 
and a free spectral range (FSR) of 20 GHz. At a red- 
shift of z = 7, a box 4 degrees on a side and 20 GHz 
deep corresponds to a rectangular box with a comoving 
angular width of 443 Mpc/h and a depth along the line 
of sight of 175 Mpc/h. However, since larger FSRs are 
easily achieved with diffraction grating architectures and 
would allow for better measurement of the reionization 
signal and separation of its foregrounds, the instrumental 
concept presented here covers the 220 GHz atmospheric 
window with an FSR of 125 GHz. Concretely, covering 
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Fig. 13. — The cross-correlation coefficient of the CII and 21-cm emission for lm and 10m aperture at z = 6, z = 7 and z = 8. The error 
bars of r are also shown (red solid), and the blue dashed ones are the contribution from the 21-cm emission with the LOFAR (left panel), 
and with the SKA (right panel). We find the 21-cm noise dominates the errors at z = 6 and 7. In the two top panels, the rcn,Hl(&) at 
z = 7 (green dashed) and z = 8 (cyan dotted) are also shown to denote to show the evolution of the ionized bubble size at these redshifts 
relative to z = 6. As expected, the bubble size is greater at z = 6 than that at z = 7 and 8. 



TABLE 1 

Experimental Parameters for a Possible CII Mapping Instrument. 



Aperture diameter (m) 


1 


3 


10 


Survey Area (As; deg 2 ) 


16 


16 


16 


Total integration time (hours) 


4000 


4000 


4000 


Free spectral range (B v \ GHz) 


185-310 


185-310 


185-310 


Freq. resolution (5 U \ GHz) 


0.4 


0.4 


0.4 


Number of bolometers 


20,000 


20,000 


20,000 


Number of spectral channels 


312 


312 


312 


Number of spatial pixels 


64 


64 


64 


Beam size a (^ eam ; FWHM, arcmin) 


4.4 


1.5 


0.4 


Beams per survey area a 


2.6 x 10 3 


2.3 x 10 4 


2.6 x 10 5 


cr p i x : Noise per detector sensitivity a (Jyy^/sr) 


2.5 x 10 6 


2.5 x 10 6 


2.5 x 10 6 


^pix : Integration time per beam a (hours) 


100 


11 


1.0 


z = 6 V pix (Mpc/h) 3 


217.1 


24.1 


2.2 


z = 7 y pix (Mpc/h) 3 


332.9 


37.0 


3.3 


z = 8 V pix (Mpc/h) 3 


481.3 


53.5 


4.8 


z = 6 P^ 11 (Jy/sr) 2 (Mpc/h) 3 
z = 7 Pg u (Jy/sr) 2 (Mpc/h) 3 


5.4xl0 9 


5.4xl0 9 


5.3xl0 9 


4.8xl0 9 


4.9xl0 9 


4.8xl0 9 


z = 8 P° n (Jy/sr) 2 (Mpc/h) 3 


4.4xl0 9 


4.4xl0 9 


4.3xl0 9 



values computed at 238 GHz, corresponding to CII at z = 7. 



from 185 to 310 GHz with a spectral resolution of 0.4 GHz 
allows measurement of CII in the range 5.1 < z < 9.3 
with a smallest redshift interval Az of 0.01. 

The integration time per beam on the sky required to 
survey a fixed area depends on a number of parameters 
including the size of the aperture, the number of inde- 
pendent spatial pixels in the spectrometer, and the band- 
width of each spectral element. Changing the survey area 
or bandwidth will affect the minimum k values probed in 
the 3-d power spectrum as well as the number of modes 
used in the error calculation. To generate concrete exam- 
ples in this discussion, we concentrate on calculating the 



CII power spectrum at z = 6, 7 and 8 and assume that 
we will make use of 20 GHz effective bandwidths at fre- 
quencies centered at these redshfits; such effective bands 
are readily achieved by summing neighbouring high res- 
olution spectral bins. At z = 7, for example, a 20 GHz 
bandwidth corresponds to a Az = 0.62. Averaging over 
a larger spectral window to reduce noise will make the 
cosmological evolution along the observational window 
non- negligible. 

To understand the effect of spatial resolution on this 
measurement we consider three possible primary aper- 
ture sizes of 1, 3 and 10 m; the resulting instrumental 
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parameters are listed in Table 1. These apertures cor- 
respond yield beams of 4.4, 1.5, and 0.4 arcmin and co- 
moving spatial resolutions of 8.3, 2.8, 0.8 Mpc/h at an 
observing frequency of 238 GHz (z = 7), respectively. 
These apertures probe the linear scales (k < 1.0 h/Mpc) 
and are well matched to 21-cm experiments. In this ex- 
perimental setup light from the aperture is coupled to a 
imaging grating dispersive element illuminating a focal 
plane of bolometers, similar to the Z-Spec spectrometer 
but with an imaging axis (Nay lor et al. 2003; Bradford 
et al. 2009). 

We assume a fixed spectral resolution of ~ 500km/s 
(= 0.4 GHz as discussed above) in each of the three cases 
giving a spatial resolution of 3.5 Mpc/h and a maximum 
k mode of k ~ 0.91 h/Mpc. Current technology can al- 
low fabrication of ~ 20, 000 bolometers in a single focal 
plane; for the instrument concept we freeze the number 
of detectors at this value. The spectrometric sensitivity 
of the grating spectrometer design is estimated using a 
model of the instrument which includes loading from the 
CMB, atmosphere and instrument contributions assum- 
ing realistic temperatures and emissivities for each. 

As an example, using equations in Appendix A, for a 3 
m aperture at 238 GHz and a NEP Bet = 1 x 10" 18 W 
HZ -1 / 2 , a value readily achieved using current detec- 
tor technology, we find a noise equivalent flux density 
(NEFD) of 0.42 Jy s 1 / 2 per pixel per spectral resolution 
element. This NEFD can then be converted to the units 
of Jy s^/sr used in Table 1 by multiplication by the 
solid angle of the telescope response. For observations 
corresponding to z = 7, this results in a noise per pixel 
of 2.5 x 10 6 Jy s^/sr. 

These noise power spectra based on Table 1 are shown 
for the 1 and 10 meter aperture in Fig. 8. The proposed 
CII experiments involve three different apertures at 1, 
3, and 10 m, but all three options make use of the same 
spectrometer with a total ~ 20,000 bolometers that make 
up 64 individual pixels on the sky. The statistical detec- 
tion shown in figure [5] using the instrumental parameters 
listed in Table 1 can be obtained by noting that the error 
on each of the binned pow er spectrum measurements is 
AP = (P™ + PcnVv^Vm, where P C n is the CII power 
spectrum, including shot-noise and N m is the number of 
/c-modes available for each power spectrum measurement. 

The noise parameters are such that with an integra- 
tion time of 4000 hours the CII power spectrum at 
z = 6 is detected with a signal-to-noise ratio, S/N = 



TABLE 2 

Sensitivities of LOFAR and SKA 21-cm experiments at 

150MHz. 



Z^bins I t 



of 12.6, 13.2, and 13.5 for 1, 3, and 



P(k) V 

AP(k)J 

10 m aperture cases, respectively. The corresponding 
values for z = 7 are 2.2, 2.2, and 2.3 respectively, while 
at z = 8 the signal-to-noise ratio values are less than 1 
for all three options. 

For the cross-correlation with 21-cm data, we assume 
that observations will be conducted in overlapping ar- 
eas on the sky by both the CII and the 21-cm exper- 
iment. This is a necessary feature of cross-correlation 
study leading to some coordination between two exper- 
iments. We reduce both data sets to a common, and a 
lower resolution cube in frequency and count the total 
number of common modes contributing to a given bin in 
k space using the sa me method for noise calculation as 
in iGong et al.l (J2011I ) . The error in a given k bin can be 



Instrument 


LOFAR 


SKA 


FoV (deg 2 ) 


25 


25 


Bandwidth (MHz) 


12 


12 


Freq. resolution (MHz) 


0.25 


0.25 


System temperature (K) 


449 


369 


maximum baseline (m) 


2000 


5000 


Total integration time (hours) 


1000 


4000 


total collecting area at 150 MHz (m 2 ) 


2.46X10 4 


1.4xl0 6 


z = 6 P^ lcm [K 2 (Mpc/h) 3 ] 


7.8X10 -1 


l.OxlO -3 


z = 7 P^ lcm [K 2 (Mpc/h) 3 ] 


1.0 


1.5XKT 3 


z = 8 P% lcm [K 2 (Mpc/h) 3 ] 


1.4 


2.3xl0- 3 



written as ^(P<5 h ,h + PciiPb) / (2N m (k)) , where N m (k) 

is the number of modes falling in that bin. 

When describing the 21-cm observations we assume 
sensitivities similar to LOFAR with parameters A tot = 
2.46 x 10 4 m 2 and T sys = 449K at 150 MHz, bandwidth of 
12 MHz, and resolution of 0.25 MHz (see Table 2 for SKA 
parameters). Here, we further assume that the baseline 
density distribution for LOFAR and SKA is constant on 
the u-v plane up to a maximum baseline D max . This is 
a reasonable approximation for the dense central core of 
these experiments which is the only configuration we're 
using for the measurements. The cross-correlation is de- 
tected with a signal-to-noise ratio of 3 to 4 at z = 6 
and 4 to 5 at z = 7 when any of the three CII op- 
tions is combined with LOFAR. There is a factor of a 
3 to 4 improvement when combined with SKA leading to 
signal-to- noise ratios of 10, 11 and 12 for the z = 6 cross- 
correlation with 1,3, and 10 m options, respectively. At 
z = 7, the signal-to-noise ratios are 10, 12, and 13, while 
at z = 8 they are 3, 4 and 5, respectively. While the CII 
power spectrum cannot be detected at z = 8 using any 
of the three options in Table 1, the cross-correlations can 
be detected when one of the proposed CII experiments 
is combined with SKA. 

The cross-correlation coefficient between CII and 21- 
cm (Fig. [T3|) shows a scale that r*cii,Hi(fc) begins to rise 
from -1 to due to the size of the ionized bubbles sur- 
rounding galaxies during reionization (Eq. [20]). This 
scale captures both the size of the typical ionized bub- 
bles and the mean value of the ionization fraction xi . The 
observational measurement of this rise-up scale is chal- 
lenging for a first-generation CII experiment involving 1 
or 3 m aperture combined with a first-generation 21-cm 
experiment like LOFAR, but is likely to be feasible with 
an improved CII experiment combined with SKA. 

8. CONCLUSION AND DISCUSSION 

In this paper, we have estimated the mean intensity 
and the spatial intensity power spectrum of the CII emis- 
sion from galaxies at z > 6. We first calculated the CII 
intensity analytically for both the ISM of galaxies and 
the diffuse IGM and showed that the CII emission from 
dense gas in the ISM of galaxies is much stronger than 
that from the IGM. Then, to check the analytical calcu- 
lation, we used the hot gas or the hot+warm gas from a 
simulation to find the CII mass in a halo, and further cal- 
culated the CII number density and intensity. We found 
that the two methods are in good agreement especially at 
high redshifts we are interested in. Next, we computed 
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the CII clustering power spectrum assuming the CII lu- 
minosity is proportional to the CII mass, Leu ~ Men- 
We compared our CII power spectrum with that derived 
from the ^cn-^co(i-o) relation, and found that they are 
consistent in the la level. We also explored the contam- 
ination of the CII emission by the CO lines at lower red- 
shift, and found the contamination can lead to 2% and 
30% enhancement of the CII power spectrum at z = 7 
and 2 = 8. 

To reduce the foreground contamination and to im- 
prove our scientific understanding of reionization we pro- 
pose here a cross-correlation study between the CII and 
21-cm emission in the overlapping redshift ranges and the 
same part of the sky. The cross-correlation exists since 
they both trace the same matter distribution. At large 
scales the correlation is due to ionized bubbles surround- 
ing CII bright galaxies while at smaller scales both CII 
galaxies and ionized bubbles trace the underlying den- 
sity field. We have outlined 3 potential CII experiments 
using 1, 3 and 10 m aperture telescope outfitted with a 
bolometer array spectrometer with 64 independent spec- 
tral pixels. A 1 or 3 meter aperture CII experiment is 



matched to a first generation 21-cm experiment such as 
LOFAR while an improved CII experiment can be opti- 
mized to match a second generation 21-cm experiment 
like SKA. We find that the overall ability to extract de- 
tails on reionization requires a careful coordination and 
an optimization of both CII and 21-cm experiments. We 
have not discussed the issues related to foregrounds, in- 
cluding AGN-dominated radio point sources that domi- 
nate the low-frequency observations and dusty galaxies 
that dominate the high-frequency CII observations and 
Galactic foregrounds such as dust and synchrotron. In 
future papers we will return to these topics and also ex- 
pand our discussion related to the instrumental concept 
as we improve the existing outline to an actual design. 
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APPENDIX 



In this Appendix we summarize a set of key equations necessary to obtain the final noise power spectrum of the 
intensity fluctuation measurements. We focus here on a single aperture scanning experiment with a bolometer array, 
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but intensity fluctuation measurements can also be pursued with interferometeric measurements (especially 21-cm 
and CO). In the context of 21-cm cross-correlation, we also provide the interferometer noise formula below without a 
detailed derivation. 

Unlike the case of C MB, where intensity fluctuation measurements are two-dimensional on the sky with the noise 
given in Knox (1995), the CII line intensity fluctuation measurements we pursue here are three-dimensional with 
information related to the spatial inhomogeneities both in angular and redshift space. For simplicity, we ignore redshift 
space distortions here and only consider a case where spatial variations are isotropic, thus the three-dimensional power 
spectrum is directional independent. The line of sight mode fen and angular/transverse mode k± are related to each 

of the k-mode via, fey = /i/c, k± = Jk 2 — k 2 , as a function of the angle in Fourier space fi. 

The noise power spectrum of intensity fluctuations take the form 

P™\k,») = V pi J^e^fH^f (1) 

where a p i x is the noise per detector sensitivity, £°^ is the integration time per beam/pixel, and Vpi x is the pixel volume 
in real space. The exponential factor captures both the spatial (A^) and radial (£4) resolution. The spatial resolution 
is set by the instrumental beam while the resolution along the radial direction are set by the spectral bin S u . These 
are respectively, £2(2) — 27r/(r(z)#beam) and k\{z) = H(z)v/[c(l + z) 5 y ] : where r(z) is the comoving radial distance 
and H{z) is the Hubble parameter. The former can be obtained from the integral r(z) = J* c/H(z) dz. For reference, 
fci ~ 0.29(5 arcmin/# beam ) (h/Mpc) and k 2 ~ 0.69(400 WRz/S u ) (h/Mpc) at z = 7. 
The total variance of the power spectrum is then 

var[P C ii(fc,M)] = [P C n(*0 + *$"(*, M)f , (2) 

where the first team denotes the usual cosmic variance. In the case where a spherically averaged power spectrum 
measurement is pursued, one can simply take the minimum variance estimate given the above /i-dependent noise 
associated with the exponential factor. Assuming a p i x is /i-independent, the /i-averaged variance of the power spectrum 
meausrement is then simply (Lidz et al. 2011) 

[Pcu(k) + P$ u (k)] 2 
var[P CII (fc)] = ^ , (3) 

where N m (k) is the number of total modes that leads to the power spectrum measurement at each k and 

PS U (k) = V pix ^| . (4) 

6 pix 

The number of modes at each k for the P(k) measurement is 

N m (k)=2nk 2 Ak-^, (5) 

where Afc is the Fourier bin size and (2tt) 3 /Vs is the resolution in Fourier space. Note that the result is obtained by 
integrated the positive angular parameter only (i.e. < /i = cos(#) < 1). 

The total survey volume is Vs = r(z) 2 yAsB u ^ where As is the survey area (in radians), B v is the total bandwidth 
of the measurement, y = A(l + z) 2 /H(z) is the factor to convert the frequency intervals to the comoving distance at 
the wavelength A. A useful scaling relation for Vs for CII measurements is 



whereA is the emission line wavelength in the rest frame. 

In Pjy H (fc), the volume surveyed by each pixel is Vpi x = r 2 yA p [ x 5 u . Here A p i x is the spatial area provided by the 
beam, or an individual pixel depending on the instrument design (in radians). A simple scaling relation for Vpi x as a 
function of the redshift z is 



2 



Y M » - 1-1 x 10 8 (Mpc/h)= (jj^-) /i±i y^-) (jjAg-) . ,7) 

Bolometer noise formula: To obtain <r p i x , for the CII case with a bolometer array, we make use of a calculation 
that estimates the noise equivalent power of the background (NEPbg) via the quadrature sum of contributions from 
Phonon, Johnson, Shot and Bose noise terms using these loadings and estimates of realistic instrument parameters 
from existing experiments (e.g. Brevik et al. 2010). Similarly, the total NEP to t of the instrument is the square root of 
the quadrature sum of NEPbg and a detector with NEPv et = 1 x 10 -18 W HZ -1 / 2 , a value readily achieved using 
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current detector technology. NEP to t does not depend on the size of the telescope aperture; to convert this to the 
corresponding noise equivalent flux density (NEFD) in Jy s 1 / 2 we use 

NEFD = Z^ Ptot (8) 

rj S y S AV2Av exp(-^ sky ) 

where A is the area of the telescope aperture. The numerical values obtained are in Table 1. 

21- cm interferometer noise formula: When estimating the cross-correlation, for simplicity, we assume that the noise 
power spectrum amplitude P^ lcm is a constant and it takes the form (e.g., Santos et al. 2010) 

«"^' = *'='4dSiM' (9) 

where Da is the comoving angular diameter distance, A to t is the collecting area for one element of the interferometer, to 
is the total integration time and the function n() captures the baseline density distribution on the plane perpendicular 
to the line of sight, assuming it i s already rotationa l ly invariant (k is the moduli of the wave mode k and is the angle 
between k and the line of sight) (Santos et al. 201 fl [Gong et al.ll2011f ). The values are LOFAR and SKA are tabulated 
in Table 2. 



